Submitted to ApJ 

Preprint typeset using lAT^^X style cmulateapj v. 11/12/01 



PROBING THE REIONIZATION HISTORY OF THE UNIVERSE USING THE COSMIC 
MICROWAVE BACKGROUND POLARIZATION 

Manoj Kaplinghat\ Mike Chu\ Zoltan Haiman^"' Gilbert P. Holder^, Lloyd Knox^ 

AND CONSTANTINOS SkORDIS^ 
Submitted to ApJ 

ABSTRACT 

The recent discovery of a Gunn-Petcrson (GP) trough in the spectrum of the rcdshift 6.28 SDSS quasar 
has raised the tantahzing possibility that we have detected the reionization of the universe. However, a 
neutral fraction (of hydrogen) as small as 0.1% is sufficient to cause the GP trough, hence its detection 
alone cannot rule out reionization at a much earlier epoch. The Cosmic Microwave Background (CMB) 
polarization anisotropy offers an alternative way to explore the dark age of the universe. We show that 
for most models constrained by the current CMB data and by the discovery of a GP trough (showing 
that reionization occurred at 2 > 6.3), MAP can detect the reionization signature in the polarization 
power spectrum. The expected l-cr error on the measurement of the electron optical depth is around 
0.03 with a weak dependence on the value of that optical depth. Such a constraint on the optical depth 
will allow MAP to achieve a 1-a error on the amplitude of the primordial power spectrum of 6%. MAP 
with two years (Planck with one year) of observation can distinguish a model with 50% (6 %) partial 
ionization between redshifts of 6.3 and 20 from a model in which hydrogen was completely neutral at 
redshifts greater than 6.3. Planck will be able to distinguish between different reionization histories even 
when they imply the same optical depth to electron scattering for the CMB photons. 

Subject headings: cosmic microwave background — cosmological parameters — cosmology: theory — 
galaxies: formation — early universe 



1. INTRODUCTION 

How and when the intergalactic medium (IGM) was 
reionized is one of the long outstanding questions in cos- 
mology, likely holding many clues about the nature of the 
first generation of light sources and the end of the cos- 
mological "Dark Age" (see Barkana & Loeb, 2001, for a 
review of our current understanding) . The lack of strong 
HI absorption (the "Gunn-Peterson " , GP, trough) in the 
spectra of high redshift quasars has revealed that the in- 
tergalactic medium (IGM) is highly ionized between red- 
shifts ^ z ^ 6. On the other hand, the lack of a strong 
damping by electron scattering of the first acoustic peak 
in the temperature anisotropy of the cosmic microwave 
background (CMB) radiation has shown that the universe 
was neutral between the redshifts 30 < z ^ 10^. Together 
these two sets of data imply that most hydrogen atoms 
in the universe were reionized during the redshift interval 
6 ^ z ^ 30. 

The recent discovery by Becker et al. (2001) of the bright 
quasar SDSS 1030-^0524 in the Sloan Digital Sky Survey 
(SDSS) at redshift z = 6.28 has, for the first time, re- 
vealed a full GP trough, i.e., a spectrum consistent with 
no flux at a substantial stretch of wavelength shortward 
of (1 -I- z)Aq, ~ 8850A. This discovery has raised the tan- 
talizing possibility that we are detecting reionization oc- 
curring near redshift z 6.3. The lack of any detectable 
flux indeed implies a strong lower limit xh ^ 0.01 on the 
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mean mass-weighted neutral fraction of the IGM at 2 ~ 6 
(Fan et al., 2002; Pentericci et al., 2002). On the other 
hand, because of the large opacities of the hydrogen Ly- 
man series, it is considerably difficult to push this method 
much further. In particular, in order to prove that we are 
probing into the neutral epoch, one would like to directly 
infer xh ~ 1- While in principle it is possible to infer this 
from the lack of any flux in a high enough S/N spectrum, 
in practice the required integration times are implausibly 
long. 

In this paper, we discuss an alternative way of probing 
deeper into the dark ages. CMB polarization anisotropy 
at large angles is very sensitive to the optical depth to elec- 
tron scattering for the CMB photons (Basko & Polnarev, 
1980; Hogan et al., 1982; Zaldarriaga, 1997; Haiman & 
Knox, 1999, and references therein). In a model with no 
reionization, the polarization signal at large angles is neg- 
ligible. However, CMB photons scattering in a reionized 
medium greatly enhance the polarization signal - making 
it very likely that such a signal (based on present CMB 
temperature anisotropy data and the GP trough) will be 
detected by the ongoing CMB satellite experiment MAP. 
We also show that Planck (and for large optical depths, 
MAP), will have the power to discriminate between differ- 
ent reionization histories even when they lead to the same 
optical depth. 

The results from CMB polarization experiments will 
deepen our understanding of the physics of reionization. 
They will lead to better constraints on the models of reion- 
ization and, thereby, tell us much more about the first 
sources of light, and indeed, the first structures to form. 

2. PROBES OF THE DARK AGE 

Very little is known about the first sources and the 
reionization of the universe. Given the existence of the 
GP trough in the spectrum of the z=6.28 quasar, it is 
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possible to argue on the basis of theoretical models and 
numerical simulations that reionization indeed occurred 
close to z = 6.3. The key ingredient in such arguments 
is that neutral fractions inferred from a sample of high 
redshift quasars show a rapid rise with redshift over the 
range 5.5 ^ z ^ 6 (Songaila & Cowie, 2002). Semi- 
analytical methods (Haiman & Loeb, 1998, for example), 
as well as numerical simulations (Gnedin & Ostriker, 1998; 
Gnedin, 1999, for example) suggest that the initial phase 
of reionization proceeds very rapidly. During the "over- 
lap" phase, when the individual HII regions of the ion- 
izing sources percolate, the mean neutral fraction drops 
from unity by several orders of magnitude in a small frac- 
tion of the Hubble time. The reason for this rapidity is 
the fast timescale for the formation rate of the ionizing 
sources, which correspond to the high-u peaks of the orig- 
inal density field. More recently, simulations of cosmolog- 
ical reionization (Ccn & McDonald, 2002; Gcndin, 2002; 
Fan ct al., 2002), geared specifically to interpret the spec- 
trum of SDSS 1030-h0524, argued that the IGM is likely 
neutral at 2; ^ 6.5. 

While there may be a theoretical bias for reionization oc- 
curring close to z = 6.3, it is conceivable for reionization to 
last over a considerably longer redshift interval. This could 
be the case if the formation rate of the ionizing sources 
does not parallel the collapse of high-cr peaks, but is more 
gradual. Complete reionization could also be delayed for 
Az > 1 due to screening by minihalos with virial temper- 
atures less than 10** K (Haiman et al., 2001; Barkana & 
Loeb, 2002). In addition, the decrease of the mean neutral 
fraction would be more gradual if the ionizing sources had 
a hard spectrum. Reionization by X -rays was considered 
recently by Oh (2001) and Venkatesan et al. (2001). In 
contrast to a picture in which discrete HII regions even- 
tually overlap, in this case, the IGM is ionized uniformly 
and gradually throughout space. 

An observational probe of the redshift history of reion- 
ization would be invaluable in constraining such scenarios, 
and to securely establish when the Dark Age ended. Al- 
though detections of the hydrogen GP trough can not pro- 
vide this type of constraint, a possibility may be to use cor- 
responding absorption troughs caused by heavy elements 
in the high redshift IGM. Recently, Oh (2002) showed that 
if the IGM is uniformly enriched by metals to a level of 
Z = 10~^ — 10~^ Zq, then absorption by resonant lines of 
01 or Sill may be detectable. The success of this method 
depends on the presence of oxygen and silicon at these 
abundance levels, and also requires that metal enrichment 
precede hydrogen reionization. 

An alternative method to probe the reionization history 
is to utilize the systematic changes in the profiles of Lyman 
alpha emission lines towards higher redshift. As argued 
in Haiman (2002), the increasing hydrogen IGM opac- 
ity towards higher redshift would make the emission lines 
appear systematically more asymmetric, and the appar- 
ent line-center systematically shifted towards longer wave- 
lengths, as absorption in the IGM becomes increasingly 
more important, and eliminates the blue side of the line. 
Because of the intrinsically noisy Ly alpha line shapes, this 
method will require a survey that delivers a large sample 
of Ly alpha emitting galaxies (Rhoads & Malhotra, 2001). 

Finally, it has been suggested that future radio tele- 



scopes could observe 21 cm emission or absorption from 
neutral hydrogen at the time of reionization (e.g., Hogan 
& Rees, 1979; Tozzi et al., 2000). This would provide a di- 
rect measure of the physical state of the neutral hydrogen 
and its evolution through the time of reionization. The 
detectability of the signal depends sensitively on the de- 
tails of reionization and will require isolation of a signal 
that is much smaller than foreground contaminants. Re- 
cently, Carilli et al. (2002) also considered the radio equiv- 
alent of the GP trough. Unlike the Lya case, the mean 
absorbtion by the neutral medium is about 1% at the red- 
shifted 21 cm. While Carilli et al. (2002) considered the 
21 cm absorbtion from large-scale structure using simu- 
lations, Furlanetto & Loeb (2002) and Ihev et al. (2002) 
have used semi-analytic methods to look at the observable 
features (in 21 cm) of minihalos and protogalactic disks. 
These studies argue that the 21 cm observations would 
yield robust information about the thermal history of col- 
lapsed structures and the ionizing background, sufficiently 
bright radio-loud quasars exist at Ze > 6.3. 

As mentioned in the introduction, an alternative method 
to probe reionization is via the CMB anisotropy. The 
advantage of this method is that it probes the presence 
of free electrons, and can therefore detect xh = 0.1 and 
Xh = 10~^ with nearly equal sensitivity. Physically, the 
CMB and the GP trough therefore probe two different 
stages of reionization: the CMB is sensitive to the initial 
phase when xh first decreases below unity, and free elec- 
trons appear, say, at redshift Zg (see Figure 2). On the 
other hand, the (hydrogen) GP trough is sensitive to the 
end phase, or, going backward in cosmic time, when neu- 
tral hydrogen atoms first appear; say, Zh- In currently 
popular theories, these two phases coincide to ^ 10 per- 
cent of the Hubble time, i.e. Ze « zh- However, as argued 
above, one can conceive alternative theories in which the 
two phases are separated by a large redshift interval, and 

Ze > ZH- 

Our goal in this paper is to quantify the power of CMB 
anisotropy to tell the difference between Ze ^ Zh and 
Ze ~ Zh ~ 6.3, and more broadly, whether there was sub- 
stantial activity in the universe prior to z=6.3. A similar 
question was addressed recently by Gnedin & Shandarin 
(2002), who consider the effects of reionization on maps of 
the temperature anisotropy. Our work is complementary, 
in that we focus on signatures of an extended period of par- 
tial reionization on the polarization anisotropy. There has 
also been some recent work on the topic of extended reion- 
ization by Bruscoli et al. (2002). Their results stressed that 
although CMB experiments will determine r, the redshift 
of reionization will remain uncertain since that depends 
on the reionization history. We point out here that in fact 
there is enough information in CMB polarization to enable 
a detailed study of the reionization history. 

In addition to probing high- redshift structure forma- 
tion, our calculations have importance for fundamental 
cosmology. In particular, the presence of free electrons 
damps the primary anisotropy, and results in a strong de- 
generacy between the amplitude of primordial fluctuations 
{A) and the electron scattering optical depth (r). The 
quantity A is one of the most important parameters in 
inflationary cosmology, as it directly probes the inflaton 
potential. We show that the knowledge we gain about r 
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from MAP data will fix A with a precision of about 5%. 

The T ~ A degeneracy and the effect of CMB polar- 
ization on it has been studied before (Zaldarriaga ct al., 
1997; Wang ct al., 1999; Eisenstein et al., 1999; Tcgmark 
ct al., 2000; Pnmct ct al., 2000; Vcnkatcsan, 2002), using 
Fisher matrix methods. The resulting measures of r and 
A have been shown to be powerful constraints on models 
of reionization (Venkatesan, 2000, 2002), and it has also 
been shown that introducing some reasonable constraints 
on models of reionization can allow improved estimates of 
cosmological parameters. In previous work, the epoch of 
reionization has been treated as a sharp transition at a 
single rcdshift, leading to all information on reionization 
from CMB measiircmcnts encoded by a single number, t. 
We show below that tlica-o is significantly more information 
in the large angle polarization anisotropy measurements. 
This will eventually allow more sophisticated tests of mod- 
els of reionization. 

3. REIONIZATION SIGNATURES IN THE CMB 

Reionization affects the CMB anisotropy in a simple 
manner. On large angular scales, reionization has no ef- 
fect, while on smaller scales it damps the anisotropy by 
a factor of exp(— 2r). The division into large and small 
scales is dictated by the angle subtended by the horizon 
at reionization. The multipole moment this angular scale 
maps onto is given by Ir = DA{z{rir))/r]r, where r]r is the 
visibility function weighted conformal time (Hu & White, 
1997) and Da{z) is the angular diameter distance to red- 
shift z. The above assumes that the photons scattering 
off of the reionized electrons do not pick up additional 
anisotropy. While not strictly correct, this approximation 
is good for r < 0.3. The electron scattering optical depth 
is given by 
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where ctt is the Thompson cross section, Ue is the elec- 
tron abundance, and uji, = Q,bh? and LOm = ^mh"^ are the 
baryon and matter densities (respectively) scaled to the 
critical density today. Eq. (1) is correct for flat w = —1 
models (Hu & White, 1997), and we have assumed com- 
plete ionization of H and "^He to H+ and ^He"*" respectively 
out to a redshift of Zg. 

The effects of reionization on CMB polarization anisotropy 
are similar to that on the temperature anisotropy. At 
I ^ Ir, the polarization angular power spectra are sup- 
pressed by exp(— 2t). However, in contrast to the tem- 
perature power spectrum, the effect at I < Ir is quite dra- 
matic (sec Figure 3). Note that there is no analog of the 
Sachs- Wolfe (or the integrated Sachs- Wolfe) effect for po- 
larization. All of the polarization signal is generated at 
the last scattering surface at z ^ 1,000 in the absence of 
reionization. Thus without reionization, the polarization 
power at small multipoles is negligible. In the reionized 
imivcrse, the rescattering of photons in the presence of a 
large quadrupole temperature anisotropy results in a broad 
peak in polarization power at low I. For a detailed discus- 
sion of this "reionization bump", we refer the reader to 
Zaldarriaga (1997). Here we provide an overview. 



We first note that since most of the optical depth to 
reionization is generated prior to the onset of curvature or 
dark energy domination, the reionization bump depends 
on curvature and dark energy only through their effect 
on Da- The polarization anisotropy is sourced by the to- 
tal quadrupole anisotropy - which for practical purposes 
is just the temperature quadrupole. The relevant quan- 
tity is the quadrupole at the reionization surface. The 
quadrupole produced at the last scattering surface (as the 
tight coupling breaks down) is small; however free-streaming 
after last scattering transfers power from the monopole to 
higher multipoles thus producing an appreciable quadrupole 
during reionization. The modes contributing to the reion- 
ization signal are A: ~ 2 /rjr which become sub-horizon well 
after recombination. Thus the quadrupole on these scales 
does not depend on the baryon density. It does depend 
on the matter density when tOm ^ 0.05 since then there's 
significant radiation around when the relevant modes en- 
ter the horizon and the integrated Sachs- Wolfe effect will 
change the quadrupole. For the most part, the reioniza- 
tion bump is only sensitive to the electron optical depth 
and the amplitude of primordial potential fluctuations. 

The arguments outlined above allow us to obtain the 
low-? E mode polarization anisotropy spectrum and its 
cross term with the temperature spectrum by simply scal- 
ing from a fiducial model. The temperature auto-correlation, 
polarization auto-correlation and the temperature-polarization 
cross-correlation spectra are denoted by C'ti , Cei and Cci 
respectively. The E subscript refers to the E (electric or 
scalar) mode and the corresponding pseudo-scalar B mode 
is absent from scalar perturbation generated polarization 
anisotropy (Kamionkowski et al., 1997; Zaldarriaga & Sel- 
jak, 1997). In terms of Cxi = 1(1 + l)Cx;/27r (where X 
stands for E,T or C), the scaling from the fiducial model 
(denoted by "*") is given by (Kaplinghat et al., 2002): 
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where Ipivot = (6000 Mpc A:pivot)/\/wm(l -I- Zrf) and the 
initial potential power spectrum is taken to be of the form 
k^P(k) = A (/c/fcpivot)""^ with fcpivot = 0.05/Mpc. The 
angular diameter distance shifting is achieved through the 
relation V = l{l'r/lr) for Cei and I' = l{l'r + 0.5)/{lr + 0.5) 
for Cci ■ 

The main features of the scaling are simple to elucidate. 
In the limit of small r, the factor of 1 — exp(— r) gives 
the fraction of scattered photons. The r^/r factor is just 
a numerical fit which makes Eq. 2 more accurate. The 
dependence on the temperature quadrupole to generate 
the polarization anisotropy signal is clearly evident in the 
n ^ 1 pivot factor. The projection formula in I reflects 
the fact that most of the action during reionization hap- 
pens close to the actual epoch of reionization (since the 
visibility function falls with increasing time). The above 
approximate scalings are good to about 10% around the 
peak of the "reionization bump" for t < 0.3. For our 
fiducial ("*") model we adopt a cosmology with r — 0.15, 

— 0.02, LOm = 0.16 and Oa = 0.65. 

Although the process of reionization is expected to be 
spatially inhomogeneous, this should have little effect on 
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the large scale polarization anisotropy discussed in the 
present study. The details of the typical patch sizes of 
reionized "bubbles" , and their spatial correlations are not 
yet well-understood theoretically (Madau et al., 1999; Cia- 
rdi et al., 2000; Gncdin, 2000; Miralda-Escudc ct al., 2000; 
Nakamoto et al., 2001). However, it appears likely that 
most of the ionizing radiation was provided by galaxies, 
(Haiman, 2001; Oh & Haiman, 2002, rather than quasars, 
for example), and the typical patch sizes are expected to be 
0(1 comoving Mpc). Inhomogcncitics can remain signifi- 
cant above the typical patch size scale because of spatial 
correlations of the ionizing sources (Knox et al., 1998; Oh, 
1999). However, these correlations fall off rapidly with 
increasing distance (approximately as l/r^'^). For the 
large angle polarization features described here the rele- 
vant modes are k ~ 2/77^; we do not expect correlations 
to significantly affect such horizon scale features (also see 
Gruzinov & Hu, 1998; Liu et al., 2001). 

4. IONIZING RADIATION AND MAXIMUM OPTICAL DEPTH 

CMB anisotropy is sensitive to both cosmological pa- 
rameters and reionization history. As a result, any prior 
knowledge about reionization can, in principle, tighten cos- 
mological constraints, and vice-versa. For the most part, 
we will adopt a conservative approach and obtain con- 
straints without assuming any such a priori knowledge. In 
this section, however, we briefly consider simplifled semi- 
analytical models of reionization, and how they can tighten 
constraints inferred from CMB data. Reionization and 
its dependence on cosmology has been studied more thor- 
oughly by Venkatesan (2002) in the context of speciflc 
reionization models. Here we restrict our analysis to ob- 
taining a theoretical upper limit to the reionization optical 
depth. In any cosmology with Gaussian seed density fluc- 
tuations, the amount of non-linear mass towards high red- 
shifts declines exponentially, and hence a relatively robust 
upper limit can be derived. 

In order to reionize the universe, the ionizing source 
must have (1) produced at least ^ 1 ionizing photon (E > 
13.6eV) per hydrogen atom, and (2) the rate of produc- 
tion of ionizing photons must be sufficient to balance re- 
combinations, and keep the hydrogen ionized. Assuming a 
uniform gas distribution at the mean IGM density, the hy- 
drogen recombination rate divided by the expansion rate 
is roughly equal to (11/(1 -|-z))^'^, which implies that re- 
combinations are inevitably significant at 2; <; 10. We 
therefore focus on the second criterion above. 

The ionizing photon production rate in the early uni- 
verse depends on several factors: the rate at which gas is 
converted into new ionizing sources; the initial mass func- 
tion (IMF) of the ionizing source population, the photon 
production rate per source (as a function of mass), and 
the fraction of ionizing photons escaping into the IGM. 
Here we make the following conservative assumptions. The 
rate of converting gas into new ionizing sources is given by 
the fraction of mass collapsed into dark matter halos with 
virial temperatures of T > 10*K (computed using the for- 
malism of Press & Schechter (1974); see discussion and 
justification in Haiman (2001)). We assume that all the 
gas in these halos turns into ionizing sources, and pro- 
duces 10^ ionizing photons per baryon. This is a conser- 
vative number; although massive 1000 Mq metal-free 



stars can produce upto ~ 100 times more ionizing photons 
(Tumlinson & Shull, 2000; Bromm et al., 2001), numerical 
simulations indicate that ^ 1% of the available baryons 
are expected to turn into ionizing sources (Abel et al., 

2000) . Wc further assume that all ionizing photons escape 
the source; in comparison, escape fractions of ~ 10% are 
measured in local starburst galaxies (Leitherer et al., 1995) 
(although larger escape fractions have been tentatively in- 
ferred for a sample of redshift z ~ 3 galaxies, Steidel et al. 
(2001)). Finally, in computing the global recombination 
rate, we assume a mean clumping factor for ionized hy- 
drogen of Cint = ('^Hii)/'^Hii = 20; this is conservative 
in comparison with numerical simulations (Gnedin, 2000; 
Gnedin & Ostriker, 1997). 

The above assumptions allow us to predict the maxi- 
mum reionization redshift, and hence a maximum optical 
depth Tmax, given a cosmology. As an example, our for- 
malism yields Zmax — 20.6 and T^ax = 0.18 for a flat cos- 
mology with ojm = 0.15, ujb = 0.02 and erg = 1. Note that 
this upper limit is significantly higher than the expected 
"ballpark" values in similar cosmologies (Haiman & Loeb, 
1998; Gnedin & Ostriker, 1998), z ~ 7 - 10. We also 
reiterate that our result for Tmax is relatively insensitive 
to our input assumptions, because of the steep decline in 
the collapsed mass fraction towards high redshifts. In our 
analysis of CMB anisotropy data in the next section, we 
will study the effects of including the constraint that the 
reionization optical depth not exceed Tmax- 

5. IS THE POLARIZATION FEATURE DETECTABLE BY 

MAP? 

In this section we examine how well MAP will be able to 
determine r. This question has already been answered for 
particular models (Eiscnstein ct al., 1999; Tcgmark ct al., 
2000; Prunet et al., 2000; Venkatesan, 2002). However, the 
answer depends on the choice of cosmological parameters, 
including t. Here we calculate the probability distribution 
for (T(r), the expected uncertainty in the measurement of 
r. assuming flat adiabatic models of structure formation 
as constrained by current CMB and quasar spectral data. 

Calculating P((T(r)) is a three-step process. Step one 
is the generation of a Markov chain of locations in the 
parameter space with the remarkable property that the 
fraction of models in a given volume of parameter space is 
equal to the posterior probability in that region; i.e., a his- 
togram of the frequency of occurence of parameter values 
in the chain ("histogramming the chain") gives the shape 
of the posterior probability distribution function (Chris- 
tensen et al., 2001). We assume a flat cosmological model 
and take our (chain) parameters to be LJb, cJd, ^A, ^, n, x 
(where r = |a:|), four calibration parameters, and a param- 
eter to account for the Boomerang beam size uncertainty 
(Netterfield et al., 2002). For this section, we assume in- 
stantaneous reionization. The data we use are from DASI 
(Halverson et al., 2002), Boomerang (Netterfield et al., 
2002), DMR (Bennett et al., 1996) and Maxima (Lee et al., 

2001) . Generation of the chain requires as many likeli- 
hood evaluations as there are chain elements (in this case 
300,000 evaluations). We use DASh (Kaplinghat et al., 

2002) which rapidly and accurately calculates C/ and per- 
forms the likelihood calculation using the offset lognormal 
approximation of Bond et al. (2000). 
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Step two is the calculation, for each chain clement, of 
"derived parameters" (parameters that are derived from 
the chain parameters). Our derived parameters include z^i, 
crE{T) and (Tc{t). The latter two are the errors expected 
on T from MAP's measurement of Cf and Cf respectively. 
We calculate z^i by inverting Eq. 1. 

We calculate (7e{t) and (Jc{t) by 
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and w = 2 X 10^^ is the expected weight-per-solid angle 
(for polarization measurements) from combining MAP's 90 
GHz and 60 GHz channels for two years. In calculating our 
results throughout the paper, we set /sky — 1 (since /sky is 
expected to be close to unity for the satellite experiments). 
Interested readers can scale our results to any value of /sky 
using Eq. 4. We calculate Cei and Cci at each point in the 
chain using Eq. 2. The partial derivatives are calculated 
by finite differencing and are accurate to ~ 10% around 
the peak of the signal. Further, the derivatives are taken 
at fixed Ar = Aexp(— 2r) since this quantity will be well- 
determined by high / CMB measurements. 

Step three is the histogramming of the chain, with addi- 
tion of various prior constraints. In this case, the priors are 
tophats implemented by simply restricting the counting of 
models to those with z^ in the desired range. 

Wc have not included effects of foreground contamina- 
tion because very little is known about polarized emis- 
sion from foregrounds in the relevant frequency range. 
Tegmark et al. (2000) have tried to model our ignorance 
of the foregrounds and use the CMB maps to learn more 
about the foregrounds. The good news from their study 
is that foregrounds can be removed, but this results in 
larger error bars for the estimated cosmological parame- 
ters. With a "conservative" estimate for the amplitudes 
of foregrounds, Tegmark et al. (2000) showed that the ex- 
pected error in r could increase by a factor of about 2 rel- 
ative to the expected error assuming no foregrounds. For 
this analysis they used all 5 MAP channels to fit the cos- 
mological and their foreground model parameters jointly. 
However, the amplitude of the polarized foregrounds is un- 
known. If this amplitude is a factor of 10 larger (i.e., the 
power is a factor of 100 larger), then Tegmark et al. (2000) 
estimate that the error on r goes up by a factor of about 
10. 

For simplicity, we assume MAP's 22, 30 and 40 GHz 
channels are used entirely for foreground removal, allow- 
ing the total weight of the 60 and 94 GHz channels to 
be used for the CMB. Also, foregrounds may significantly 
contaminate Cci but not Cei, ot vice versa and hence the 
existence of two different statistical channels to determine 
T increases our chances. Note that the errors asir) and 



cc{t) are correlated and this must be taken into account 
if one were to combine them into a final error for r. 

We show the results of the exercise detailed in this sec- 
tion in Fig. 1. The four panels in the figure show contours 
of equal probability (1 and 2-a) in the (JEiT) and i7c{t) 
vs. r and Zg planes. We have plotted these contours for 
different priors on Zg- If reionization is indeed instanta- 
neous, then the 6 < Ze < 7 prior (solid contours) is the 
interesting case to study. However, as we stressed before, 
partial reionization is possible at redshifts z,. > 6, and 
hence we have included contours with more conservative 
priors: 6 < < 10 (dashed contours) and 6 < Zg (dotted 
contours). 
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Fig. 1. — 1 and 2-a joint probability contours of optical depth r 
(or Ze) and the expected uncertainty with which MAP can measure 
T in two years of operation. Measurements of r using the E (EE) 
and the C (TE) channels arc considered separately. The left panels 
show the T, crE,c{'^) contours and the right panels show the Ze, 
cij) contours. Wc have included three different priors on z^- 
The dotted contours arc for Ze > 6, the dashed contours are for 
6 < 2e < 10, and the solid contours are for 6 < Ze < 7. 



The result to take away from Figure 1 is that MAP will 
be able to determine r to about 0.02 to 0.03 regardless 
of the value of r. The lack of dependence of (Tb(t) on r 
is due to two reasons. First, at its peak Cei « r (that 
this is approximately true can be ascertained from Eq. 2) 
and hence around peak power dCEi/dr is approximately 
constant. Second, ACei is dominated by the noise 
term rather than the Cei sample variance term. The cross- 
correlation power spectrum Cci varies more slowly (than 
linear) with r and hence (Tc(t) increases with r. 

Figure 1 also includes the constraint that the optical 
depth not exceed Tmax calculated as outlined in Section 4. 
The effect of this Tmax cut is to exclude some of the large 
optical depth, r ^ 0.2, models. Without the Tmax cut, 
the l-cr contours extend to about t — 0.3 for (7e{t) (and 
T = 0.25 for acir)). With the cut, both the l-cr contours 
are restricted to t < 0.2. As mentioned in Section 4 the 
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parameters we use in our semi-analytical model to calcu- 
late Tmax are very conservative. If the optical depth is in 
fact large, then reionization models will provide significant 
complementary constraints. 

One can look at the dctcctability of the signal in another 
way. Let's denote the value of the maximum Ce,i for / be- 
tween 2 and 20 by C;,peak- One can then include C;,peak 
as a chain parameter and obtain its probability distribu- 
tion. A rough estimate of the detectability of the signal 
can be obtained from looking at if C;.poak (signal) is larger 
than the detector noise contribution; for MAP this noise 
contribution is •w~^ = 5 x 10~^^. 

The distribution of C;,pcak peaks at 5.7 x 10"-^^. 90% 
of the models with z,. > 6.3 in the chain have Q^peak > 
5.5 X 10~^^, which tolls us that its likely MAP will detect 
the reionization feature. If we now apply the Tmax cut, 
then C,,peak > 4.9 x lO'^^ for 90% of the models. If we 
further restrict Ze to between 6.3 and 7, then 90% of the 
models have Cj^peak > 3.2 x 10~^^. 

We end this section by noting that since will be 
determined with very high precision by MAP, the uncer- 
tainty in A will be dominated by the uncertainty in r. If 
MAP detects the reionization feature as discussed above, 
then the uncertainty in A is expected to be o-(ln^) = 
2cr(T) = 0.04 to 0.06. Note that this is far better than what 
can be achieved using large angle temperature anisotropy 
(COBE) data. 
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Fig. 2. — The evolution of the ionization fraction as a function 
of redshift (z) in a two-step reionization modeh ne and nji are the 
number density of free electrons and hydrogen nuclei, respectively. 
The optical depth to the last scattering surface in this model is 
0.05 and the epoch when electrons first appear is chara<;terized by 
Ze = 15 and Xe = 0.23. 



6. IS A GRADUAL TRANSITION DISTINGUISHABLE FROM 
INSTANTANEOUS? 

The answer to this question depends on the details of 
how reionization took place. For simplicity, here we as- 
sume a two-step reionization process (see Figure 2 for the 
evolution of the ionization fraction in one such model). 
We assume that at some Ze > zn — 6.3, the universe 
was partially ionized, the ionized fraction being Xe < 1.08. 
The ionized fraction is defined here as the ratio of number 
density of free electrons to that of hydrogen nulei ®. At 
around z = 6.3, new sources turned on, leading to almost 
complete ionization. If reionization happened gradually 
(and not in our idealized two-step fashion) then Xg pro- 
vides us with an estimate of the average ionized fraction 
before z = 6.3, and Zg an estimate of when the first ioniz- 
ing sources turned on. As discussed above, the large scale 
polarization anisotropy should be insensitive to the topol- 
ogy of the partially ionized regions. 

The power spectra used in the analyses described in this 
section were calculated using the publicly available CMB- 
fast code (Seljak & Zaldarriaga, 1996) which we modified 
to take into account the two-step reionization process. 

We push the arguments outlined in Section 4, to its 
limit, and assume a minimum virial temperature of 100 
K for halos to cool (via H2 molecules; see Haiman et al., 
2000) and form astrophysical objects that can contribute 
to reionization. Note that gas colder than 100 K can not 
form H2 and cool, and therefore can not contract or frag- 
ment in the dark halo. We find that this cooling cutoff 
leads to the result: Zg ^ 40. In order to have Zg exceed 
this value, the presence of non-gaussianity in the primor- 
dial power spectrum would be required. 

®The present ionized fraction is slightly larger than unity because 
of the contribution from the complete reionization of ^He — ♦ '*He"'". 
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Fig. 3. — Effect of different reionization histories on t he CMB 
anisotropy. Panel (a) shows Cei, and (b) shows r; = Cci/^/CeiCti- 
Both panels have the same set of reionization histories. The solid 
curve shows the result of no reionization. The other 3 curves have 
zh = 6.3. The dashed curve has Xe = which implies r = 0.033 
(note that is the ionized fraction for z > zjj). The dotted curve 
has r = 0.1 and Ze = 30 which requires x^ = 0.26. The dot-dashed 
curve has T = 0.1 and z^ = 15 which requires Xe = 0.89. We have 
added cosmic variance error bars to the z^ = 15 Cei (dot-dashed) 
curve. 
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Figure 3 shows the effeet of different reionization his- 
tories. The polarization "bump" due to reionization is 
unmistakable, as is the lack of power in the model with 
no reionization (solid curve). The exp(— 2t) suppression 
at large / is not very evident because of the logarith- 
mic scale for the E mode polarization spectrum in panel 
(a). Note that the variable related to cross-polarization 
plotte d in pan el (b) is the coefficient of correlation r; = 
CciNCeiCti (de Oliveira-Costa et al., 2002). 

The curves to compare in Figure 3 arc the dot-dashed 
and dotted ones. Both the curves have r = 0.1; the dotted 
curve has Ze = 15 while the dot-dashed curve has Ze = 30. 
The main effect of increasing Zg at fixed t is to spread 
the total reionization induced power over a larger range 
of /. In fact, when Zg and zh are well separated, two 
reionization "bumps" can be discerned. These are directly 
related to the bimodal nature of the visibility function for 
our 2-stcp rcionzation process. The change in the shape of 
the polarization power spectra due to partial rcionzation 
at z > Zh allows two models with the same optical depth 
but different reionization histories to be distinguished from 
each other. For reference we have plotted the cosmic vari- 
ance error bars given by ^/2/{2l + 1)Ce,i for the Ze = 15 
Ce,i curve. 

An interesting question, which is part of our broader 

goal of distinguishing between models with differing ion- 
ized fractions beyond z — 6.3, is whether one can distin- 
guish a completely neutral hydrogen content at z > 6.3 
from a partially ionized one. We outline the method used 
to compare models with different reionization histories be- 
low. 

Given the spherical harmonic coefficients derived 
from data, one can write down the likelihood function for 
a particular model (Cxi) given the data as: 



-21n(£)=^(2/ + l)/,ky 



Tl 



El 



where C^; = J2r^ 



[C'ei 

"Zml 



W 

/{2l 



] Cti 
1). Eq 



(5) 



5 was derived 



under the assumption that arc gaussian distributed 

with mean zero. We have set the damping due to the fi- 
nite experimental beam width to be unity (since we are 
only interested in the low I region), and omitted the noise 
term for Cti (since it's negligible in comparison to Cti)- 
We have no data in hand. Instead, we will assume that 
C^j is given by some fiducial model, C^i (in other words 
we perform an ensemble average over ln(£)), such that 



r<fid 



-w^ with = 0. We then look at the like- 



lihood of other models in this parameter space. Note that 
all the discriminatory power comes from I < 30. In this 
regime (as we have shown in Eq. 2) the only cosmological 
parameters the spectra are sensitive to are the amplitude 
of power spectrum, tilt and the optical depth. 

To fix the amplitude we again use the fact that the CMB 
constrains the combination At = Aexp(— 2t) very well - 
to about 10% (at l-tr) with present data (without including 
massive neutrinos). This will greatly improve with MAP. 
Note that this particular combination of the amplitude of 
the power spectrum and optical depth is not affected by 
the tensor contribution to the anisotropy since the con- 
straints come from the small angular scales. The error in 
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Fig. 4. — Constant likelihood contours at 10% (thick curves) 
and 1% (thin curves) of the maximum likelihood. The maximum 
likelihood occurs at the fiducial model, which for this plot is one 
with no ionized hydrogen at z > 6.3. The upper panel assumes two 
MAP and Planck channels as described in the text. The lower panel 
assumes one channel per experiment. The solid curves are for MAP. 
The dashed curves arc for Planck. The dot-dashed curves arc for a 
hypothetical SuperPlanck (described in the text) . The dotted curve 
is a constant t = 0.06 contour. 



At when propagated to the Cis is much smaller than cos- 
mic variance on the scales that we consider {I < 30). This 

is simple to see: since the Cis scale linearly with A, for 
any given r the fractional theoretical error in Cxi due to 
the uncertainty in At is a (At) /At. Thus even with cur- 
rent data this fractional error becomes larger than cosmic 
variance only at I > 100 - safely outside the regime we are 
interested in. 

We set the tilt to unity. Future CMB data sets will de- 
termine the tilt well - the expected errors are a{n) = 0.05 
for MAP and (T{n) = 0.02 for Planck (allowing for small 
variations of the spectral index with scale) (Eisenstein 
et al., 1999). Relaxing this constraint will not affect our 
results. For Cei and Cci, the predominant effect of tilt 
(at small I) is just a rescaling of the spectra - degenerate 
with the amplitude of the power spectrum. The fractional 
change in Cc,ei due to a non-zero n — 1 is ~ 4.5(n — 1). 
Given the Z-range of the signal and the precision with 
which n will be determined, this is not an issue. 

The experiments we consider here are MAP and Planck. 
To be conservative, as mentioned, we only assume a max- 
imum of two MAP or Planck channels. The best guess at 
the two relatively foreground-free Planck channels are the 
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Fig. 5. — Constant likelihood contours at 10% (thick curves) and 
1% (thin curves) of the maximum likelihood - which occurs at the 
fiducial models labeled by A, B and C and denoted by asterisks. 
A has Ze = 32 and Xe = 0.1; B has = 14 and Xe = 0.4; C 
has = 26 and x'e = 1.0. All the panels have the same fiducial 
models. The solid curves correspond to model A, the dashed curves 
to model B and the dot-dashed curves to model C. The upper panel 
is for a cosmic variance limited experiment, "SupcrPlanck" (see text 
for details). The middle panel assumes 1 Planck channel sensitivity 
while the lower panel is for two channel MAP sensitivity. The dotted 
curves in the lower panel are curves of constant t for r = 0.063 and 
T = 0.25. 



143 and 100 GHz channels (Tegmark et al., 2000). The 
weight-per-soUd angle (for polarization measurements) of 
the 143 GHz channel for one year of observation is ex- 
pected to be 1.35 X IQ-'^^ and from combining the 100 and 
143 GHz channels it's expected that w = 1.67 x 10^*^. 

We will not consider here the balloon borne and terres- 
trial polarization experiments since these experiments will 
not have the required sky coverage to obtain spectra at 
the large angular scales at which reionization imprints its 
signal. 

Figure 4 shows constant likelihood contours for the case 
where the fiducial model has only neutral hydrogen at z > 
6.3. The constant likelihood contours are drawn at 10% 
and 1% of the maximum likelihood. We have considered 
two separate cases plotted in the upper and lower panels 
corresponding to sensitivities equal to that of MAP and 
Planck assuming two channels or a single channel. This is 
important for MAP since each channel has equal weight. 

In the lower panel of Figure 4 we also plot the results 
for an experiment we dub "SuperPlanck" with a sensitivity 



which is 25 times better than Planck, i.e., w = 3.4 x 10""^^. 
For reference, the raw sensitivity of Planck is a 

factor of 11.6 better than MAP. Clearly, the results of us- 
ing SuperPlanck are remarkable. We have chosen Super- 
Planck to be 25 times more sensitive than Planck because 
that makes SuperPlanck nearly cosmic variance limited 
for Ze < 30. The sensitivity of SuperPlanck also makes it 
a great tool to hunt for primordial (tensor) B-mode po- 
larization signal. This signal would be detectable by Su- 
perPlanck if the ratio of tensor to scalar pertuburbations 
is larger than 0.001 (Knox & Song, 2002; Kesden et al., 
2002) which implies an energy scale of inflation greater 
than 6 x lO^^GeV. We study the results of Figure 4 in 
more detail below. 

Let us first concentrate on MAP. Using the 10% likeli- 
hood contours, we see that MAP should be able to dis- 
tinguish a completely neutral hydrogen content aX z > 6.3 
from a case where hydrogen is 50% ionized between the 
redshifts of zh = 6.3 and = 20. We see that in this 
region of parameter space the contours are essentially par- 
allel to the constant r contours - implying there is very 
little information except for the value of r itself. 

Planck, of course, can do a much better job. Planck can 
distinguish between 10% ionized fraction between zh and 
Ze = 15 and the fiducial model which has a;e = at z > 
zh = 6.3. The hypothetical SuperPlanck could distinguish 
the fiducial abrupt-reionization model from one in which 
Xe = 0.05 between zh and z,, = 15. 

We note that if the universe was indeed reionized close 
to z = 6.3, then MAP will be able to show that Ze < 13, 15 
at 2, 3-(T. Clearly, this will not be sufficient to show that 
the universe was really abruptly reionized close to z = 
6.3. However, for the same assumptions, the upper bound 
from Planck will be Zg < 7.5, 8. This will be a strong 
indication that most of the reionization took place close 
to z 6.3, although delays of Az ~ 1 will still remain 
a possibility. These results can be read oif from Figure 4 
(the 10% and 1% contours projected to the Xe = l line 
aproximately give the 2 and 3 a bounds), or derived from 
Figure 1. From Figure 1 we see that it(t) ~ 0.025 (for 
MAP) which implies (T(ze) ~ (2/3)ze(T(T)/T = 3 for Zg = 
6.3. If we assume that both MAP and Planck are noise 
(and not cosmic variance) limited for Zg = 6.3, then we 
have a{ze) — 0.4 for Planck. This is somewhat better 
than what Figure 4 implies because cosmic variance is not 
completely negligible. 

One could also ask the following question: if the real 
universe was completely reionized at some z > 6.3, then 
at what z does it become possible to be certain (say at 3-cr) 
that Ze > 6.3? We can use the discussion of the previous 
paragraph to answer this question. Taking r ~ 0.002z''^'^, 
we have that a{ze) — 333c7(t)/^/z^ and to a good approx- 
imation we can assume that (j(r) is independent of r for 
the optical depths we are considering here. Thus if the 
universe was reionized at a redshift of 13 (8), then MAP 
(Planck) will be able to claim that Zg > 6.3 at the 3-cr 
level. 

We take the above discussion to its natural conclusion 
by answering the following question: what information can 
one obtain if the universe is partially ionized at z > 6.3? 
Figure 5 is similar to Figure 4 except we now plot likeli- 
hood contours for different fiducial models. Specifically we 



9 



consider 3 fiducial models: (A) 2e = 32 and Xe = 0.1, (B) 
Ze = 14 and Xe = 0.4, (C) Ze = 26 and Xe = 1.0. Models 
(A) and (B) have the same optical depth of t = 0.065, 
while model (C) has a large optical depth of r = 0.25. We 
see that at the larger optical depth, MAP is sensitive to 
more than just r. The distribution of power in I which is 
decided by Ze is sufficiently different that the contour does 
not follow the constant r contour along the entire parame- 
ter space. For Planck this happens at the much lower r of 
0.065. It is clear from the Planck contours (upper panel) 
that cases (A) and (B) can be distinguished despite their 
equal optical depths. 

This is an interesting result. It opens up the possibil- 
ity that we can observationally determine the reionization 
history of the universe. The usefulness of suc;li data will 
be directly determined by the sophistication of the physi- 
cal modeling of the end dark ages is. We end this section 
with the tantalizing suggestion that with SuperPlanck, it 
should be possible to map the evolution of the ionized hy- 
drogen fraction Xe{z) by breaking up the za-Ze redshift 
interval into several bins and tracing the evolution of Xg- 

7. CONCLUSIONS 

The remarkable discovery of Gunn-Peterson trough in 
the spectrum of the redshift z — 6.28 quasar has brought 
into sharp focus the questions of how and when the uni- 
verse was reionized. While representing a significant new 
discovery, the power of the hydrogen Gunn-Peterson trough 
to probe the dark age is limited, because a small neutral 
hydrogen fraction (xh ~ 10^"^) is enough to give rise to 
such a trough. GMB polarization is the most promising 
way to probe deeper into the dark age; our goal in this 
paper has been to quantify this statement for forthcoming 
CMB anisotropy experiments. 

Contaminating astrophysical (foreground) sources of po- 
larization may be present at significant levels, but there is 
very little we know about them. In deriving our results, 
we have used a maximum of two channels per experiment, 
with the remainder assumed to be used as foreground mon- 
itors. Foregrounds might be more detrimental than we 
have assumed; the only way to find out is to do multi- 
frequency CMB experiments. 

Using current data we have shown that MAP will most 
likely detect the reionization feature in the polarization 
data. We have shown that MAP can measure the electron 
optical depth to 0.02-0.03 (1-cr), regardless of the value of 
the optical depth and that this measurement is sufficient 
to break the degeneracy between the optical depth and 
the amplitude of primordial density perturbations. Fur- 
ther, MAP (with two years of observation) will be able to 
distinguish a model with only neutral hydrogen at z > 6.3 
from one that was about 50% ionized between z = 6.3 
and 20. Planck, with its much higher projected sensitiv- 
ity (and one year of observation), will be able to distin- 
guish an ionized hydrogen fraction of about 6% between 
the redshifts of 6.3 and 20 from the model with no ionized 
hydrogen at z > 6.3. More strikingly, Planck (and MAP if 
the optical depth is larger than about 0.2) should be able 
to differentiate among reionization histories which lead to 
the same optical depth for CMB photons. These results 
will enormously deepen our understanding of the physics 
of reionization. 
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